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Abstract 

Many theories beyond the Standard Model predict the existence of colored scalar states, known as sgluons, lying in 
the adjoint representation of the QCD gauge group. In scenarios where they are top-philic, sgluons are expected to be 
copiously pair-produced at the LHC via strong interactions with decays into pairs of top quarks or gluons. Consequently, 
sgluons can be sought in multijet and multitop events at the LHC. We revisit two LHC Run I analyses in which events 
featuring either the same-sign dileptonic decay of a four-top-quark system or its single leptonic decay are probed. 
Adopting a simplified model approach, we show how this reinterpretation allows us to extract simultaneous bounds on 
the sgluon mass and couplings. 
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1. Introduction 

Despite its success in describing all experimental high- 
energy physics data, the Standard Model (SM) of par¬ 
ticle physics leaves many important and conceptual is¬ 
sues unanswered. As a consequence, many theoretical 
frameworks extending it have been developed over the last 
decades, and new phenomena have been searched for ex¬ 
perimentally. Weak scale supersymmetry, and in particu¬ 
lar its minimal realization known as the Minimal Super- 
symmetric Standard Model (MSSM) [Tl[2], is one of the 
most studied of those beyond the SM setups. It is however 
more and more constrained by data, and especially by the 
recent results of the LHC experiments uni- There are 
nevertheless large varieties of alternative non-minimal su¬ 
persymmetric models that deserve to be investigated and 
whose signatures may be different from the expectations 
of the minimal choice. 

Along these lines, we focus on N = 1/N = 2 hybrid [S1[S] 
and R-symmetric [7HS1 supersymmetric theories that both 
predict extra scalar partners to the SM gauge bosons. 
These additional degrees of freedom lie in the adjoint rep¬ 
resentation of the gauge group and are indeed not present 
in the MSSM. Among the new particles, the colored states 
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commonly dubbed sgluons have received special attention 
as they are expected to be copiously produced at hadron 
colliders [TUHIHl- Those fields however appear not only 
in supersymmetry but also in vector-like confining theo¬ 
ries [TS] and extra-dimensional models [501. Motivated by 
typical sgluon signatures that are similar in all these mod¬ 
els, we adopt a simplified model approach describing the 
dynamics of a scalar field lying in the octet representation 
of SU{3)c and interacting with the SM [T3|2T]. This sub¬ 
sequently allows both for a model-independent approach 
and a simplification of the non-minimal supersymmetric 
parameter spaces that in general contain hundreds of free 
parameters. 

This simplified approach has been used experimentally 
in order to search for hints of sgluons in LHC collision 
data at center-of-mass energies of 7 TeV and 8 TeV. As 
a consequence of null results, limits on the sgluon mass, 
ms, have been extracted after probing the production of a 
sgluon pair that decays into either a four top-quark system 
yielding a same-sign dilepton final state [22], or a four- 
jet final state [53]. In the former case, ms is bound to 
be larger than 800 GeV when we assume that the sgluon 
always decays into a top-antitop system. In the latter 
case, limits turn out to be weaker (ms > 300 GeV) and 
are obtained after assuming that the sgluon always de¬ 
cays into a dijet state. Additionally, stronger constraints 
could be derived in the context of single sgluon produc¬ 
tion, the sgluon mass being pushed in this case above 
the multi-TeV scale [53] |5S] . These bounds are however 
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very model-dependent and could be evaded as soon as the 
sgluon is allowed to couple to the top quark, as in realistic 
ultraviolet-complete setups such as those mentioned. We 
therefore consider a framework where sgluons can couple 
to both quarks and gluons, and then revisit constraints 
on sgluon simplified models by reinterpreting recent LHC 
analyses of all data recorded at a collision center-of-mass 
energy of 8 TeV. More precisely, we consider two CMS 
studies of four-top-quark topologies, a first one focusing 
on same-sign dilepton events |26j and a second one on sin¬ 
gle lepton events [27]. We hence derive, for the first time, 
limits on the sgluon mass and coupling strengths to the 
SM particles simultaneously. 

The rest of this paper is organized as follows. In Sec¬ 
tion]^ we briefly describe our simplified theoretical frame¬ 
work for sgluon phenomenology at the LHC, and present 
the sgluon mass dependence of both the total sgluon-pair 
production rate and the sgluon branching ratios. The rein¬ 
terpretation of the LHC analyses of Refs. El HZ] is de¬ 
tailed in Section and our conclusions are presented in 
Section |4j 

2. A simplified model for top-philic sgluon phe¬ 
nomenology 

In our study of sgluon production and decay at the 
LHC, we rely on a minimal extension of the Standard 
Model allowing for a general description of sgluon dynam¬ 
ics. To this aim, we construct a simplified model in which 
we supplement the Standard Model by a single real color- 
octet scalar held 5'“, the superscript a indicating an adjoint 
color index. The kinetic and mass terms associated with 
this held are given by the Lagrangian 

C = ^D^S^D>^Sa - ^mlS^Sa , (1) 

which includes the gauge interactions of a sgluon pair to 
gluons through the QCD-covariant derivative, 

= . ( 2 ) 

In our conventions, gg denotes the strong coupling con¬ 
stant, fbc°‘ the structure constants of the SU{3)c gauge 
group and represents the gluon held. 

Furthermore, in order to allow the sgluon for singly 
coupling to the Standard Model degrees of freedom, we 
introduce the effective Lagrangian 

Geff = + + + 

(3) 

Its hrst term consists of dimension-four interactions of the 
sgluon with a pair of top-antitop quarks whose left-handed 
and right-handed coupling strengths are denoted by 
and af, respectively. The second term of £eff models 
single sgluon interactions to gluons through a dimension- 
hve operator, the dimensionless coupling strength Ug be¬ 
ing suppressed by the theory cutoff energy scale A. In 
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Figure 1: Total cross section for sgluon pair production at the LHC, 
for a center-of-mass energy of 8 TeV, given as a function of the sgluon 
mass. The lower panel shows the sgluon decay branching ratio into 
a top-antitop pair for different coupling strengths at. 

our notations, the Ta matrices stand for the generators of 
SU{3)c in the fundamental representation, while da^'^ are 
the symmetric structure constants of the group and Pl/r 
denote the chirality projection operators. Although sglu¬ 
ons can in principle also couple in a flavor-changing-neutral 
way to different quark species, we only retain their flavor- 
conserving interactions with top quarks. While not gen¬ 
eral, this choice is motivated by minimally flavor-violating 
i?-symmetric supersymmetric models where single sgluon 
interactions are loop-induced by squarks and gluino in a 
way such that only interactions involving a pair of gluons 
or top quarks are non-negligible m- 

This setup corresponds to the scenarios of class H in¬ 
troduced in Ref. in which the sgluon is top-philic and 
only allowed to decay into a top-antitop pair or into two 
gluon-induced jets. We define our reference benchmark 
scenario by 

af = of = at = 1.5 X 10"^ , ag = 1.5 x 10"^ 

A = 1 TeV , (4) 

and impose the sgluon mass ms to lie in the [350, 900] GeV 
window. For the phenomenological investigations performed 
in this work, we additionally study deviations from this 
reference scenario by varying the o* and ag parameters in 
the ranges [0.5,5] x 10“^ and [1.35,1.65] x 10“^, respec¬ 
tively. As shown in Ref. m, those numerical values can 
be obtained in ultraviolet-complete supersymmetric mod- 
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els featuring colored superpartner masses of about 1 or 
2 TeV. 

The scenarios under consideration exhibit an enhance¬ 
ment of the production rate, at the LHC, of events con¬ 
taining four top quarks. This is illustrated on Figure [T] 
for proton-proton collisions at a center-of-mass energy of 
8 TeV. First, we present, on the upper panel of the fig¬ 
ure, the TOg-dependence of the sgluon-pair total produc¬ 
tion cross section evaluated at the next-to-leading order 
(NLO) accuracy in QCD [161 [H]. Second, we show, in the 
lower panel of the figure, the ms-dependence of the sgluon 
branching ratios into a top-antitop system for several val¬ 
ues of at- In order to calculate the sgluon-pair total cross 
section at NLO in QCD, we have followed the procedure 
described in Ref. [18] . More specifically, we have employed 
the FeynRules package [28] and its NloCT module [29| 
to generate a UFO library [30] suitable to be used within 
the MADGRAPH5_aMC@NLO framework [3T]. The cen¬ 
tral curve on Figure[2is then obtained by fixing the renor¬ 
malization and factorization scales to ms and using the 
NNPDF 2.3 set of parton distributions [32], while the un¬ 
certainty band has been derived by varying the two un¬ 
physical scales by a factor of two up and down with respect 
to ms (the inner range) and by using all parton density 
replicas provided by the NNPDF Collaboration (the outer 
range). 


3. Constraining top-philic sgluons from LHC Run 

I data 

As sketched on Figure the production of a pair of 
top-philic sgluons, whose dynamics are described by the 
model presented in Section]^ leads to an enhancement of 
LHC events containing either two top-antitop systems, or 
two dijet systems, or one of each. Since in our class of real¬ 
istic scenarios, the sgluon couples to both top quarks and 
gluons, the existing limits on its mass [22] [24] [25] do not 
directly apply and must be carefully reinterpreted. The 
most stringent and model-independent constraints have 
been derived in the context of an ATLAS analysis of par¬ 
tial LHC Run I data [32] . Sgluons are in this case searched 
for in same-sign dilepton events arising from the decay of 
a four-top (titi) system, since such a final state has the 
advantage to allow for a good experimental precision due 
to very low SM background. Instead of recasting this AT¬ 


LAS analysis, we focus in Section 3.1 on its CMS counter¬ 
part [3^ which benefits from the entire collision dataset at 
a center-of-mass energy of 8 TeV. Additionally, we choose 
to also explore the tttt topology via its single lepton plus 
jets decay channel. Although this mode exhibits consider¬ 
able but well-understood background, it allows us to ex¬ 
ploit a larger signal branching fraction. We reinterpret in 
Section |3.2| the SM tttt single lepton plus jets analysis of 
CMS [27] that also covers all 8 TeV data. Our results, that 
consist of the first attempt to reinterpret LHC analyses of 
the full 8 TeV dataset in the aim of constraining realistic 


sgluon models, are presented in Section 3.3 


3.1. Dilepton analysis 

The first topology chosen to be investigated in this 
work consists of a same-sign dilepton signature, for which 
CMS has analyzed the entire recorded dataset of collisions 
at a center-of-mass energy of ^/s = 8 TeV [33]. This CMS 
analysis consists of a collection of counting experiments 
searching for new physics from events containing two iso¬ 
lated same-sign charged leptons (ee, e/i, p,p) and jets. Re¬ 
sults are presented in multiple search regions which are 
determined by requirements on the number of selected jets 
(A(jets) with transverse momentum pT > 40 GeV and pseu¬ 
dorapidity satisfying \r]\ < 2.4, as well as on the number of 
jets identified as originating from the fragmentation of a b- 
quark (A^b-jets)) the scalar sum of the px of the selected jets 
(iLx) and the missing transverse energy {^t)- Without a 
knowledge of the correlations between the uncertainties 
on the background expectations in the search regions, the 
reinterpretation of a combination of search regions is not 
possible. Therefore we utilise the results from the search 
region 28 (SR28) only. This search region is chosen as its 
requirements can be fully emulated using the selection ef¬ 
ficiency parameterizations included in the publication and 
closely correspond to the four-top-quark signature. This 
region is defined by requirements of A^jets > 4, A^b-jets > 2, 
iLx > 400 GeV and > 120 GeV, imposed together with 
the demand of two same-sign leptons with px > 20 GeV 
and IpI < 2.4. Events containing an opposite-sign same- 
flavor lepton pair with an invariant mass Mu satisfying 
either Mu < 12 GeV or 76 GeV < Mu < 106 GeV are 
also vetoed. This was implemented in the original analy¬ 
sis in Ref. |26j to suppress background events arising from 
the decay of a low-mass bound state or multiboson pro¬ 
duction. 

The analysis of Ref. |26j contains very useful informa¬ 
tion allowing for reinterpretation studies, following closely 
the recommendations of Refs. [33] [M]. In particular it in¬ 
cludes parameterizations of the selection efficiencies which 
allows the signal acceptance to be estimated from generator- 
level information without the need for full detector simu¬ 
lation. Using these parameterizations we define an event- 
by-event weight 


^event — ^ ^ ^>2b-tags ^ ^>1SS2L , (5) 

that uses the individual efficiencies for the hadronic activ¬ 
ity Ht {sht)^ the missing energy (e^^), the selection 
of at least two 5-jets (e>2b-tags) and at least one same-sign 
lepton pair (e>iss2L), this last efficiency being obtained af¬ 
ter considering all possible permutations among the event 
leptons. In addition, a mistagging rate of a jet originating 
from a light quark or a gluon as a 5-jet of 1% is incorpo¬ 
rated in the calculation of e>2b-tags- 

Using the MADGRAPH5_aMG@NLO generator [3T| . 
we simulate the production, at the LHC with y/s = 8 TeV, 
of a sgluon pair that decays into a tttt final state. Our 
event sample is generated inclusively in the top decays, and 
we subsequently match the hard-scattering events to the 
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Pythia 6 parton showering and hadronization [3S]. Jets 
are reconstructed by using the anti-ZcT algorithm with a ra¬ 
dius parameter set to i? = 0.5 [35] as included in the Fast- 
Jet package m, while the event selection and reweighting 
procedure are implemented within the MadAnalysis 5 
framework [SHIES]- A signal acceptance of 0.60% has been 
calculated for SM tttt events using the efficiency modeling 
of Eq. ([^. This can be compared with the acceptance of 
0.49% obtained by CMS for the same process (23 . These 
two efficiencies agree within the 30% uncertainty quoted 
on the efficiency model. 

Scanning over the parameter space introduced in Sec¬ 
tion]^ we compute signal event yields for different {ms, at) 
values by using the relevant NLO cross section (see Fig- 
ure[^, a luminosity of 19.6 fb“^ and the signal acceptance 
derived with the efficiency model. The obtained results are 
then compared to the 2 (2.21) events observed (expected) 
by CMS in the SR28 region. Considering a 30% uncer¬ 
tainty on the signal yield (see above), we use the asymp¬ 
totic CLs method uni as implemented in the RooStat 
package |4T| to calculate a 95% CL observed (expected) 
upper limit on the number of signal events N in SR28. 
We have computed upper limits of 4.68 (4.89), therefore 
we exclude the region of the {ms, at) plane where N is 
predicted to be larger than 4.68 events. 

3.2. Single lepton analysis 

The CMS titi single lepton analysis has been performed 
assuming the specific kinematics of the SM tttt production, 
and it places, at the 95% CL an upper limit of 32 fb on 
the associated cross section EH- For certain values of ms 
and at, sgluon-induced contributions could enhance the 
tttt total cross section to a value larger than 32 fb, so 
that constraints on the parameter space could be extracted 
from the tttt single lepton analysis of CMS. However, in 
order to properly estimate which regions of the parameter 
space are disfavored by data, it is necessary to verify that 
the final state kinematics in the decay of a sgluon pair are 
similar to the SM case. To this aim, we employ a simplified 
matrix-element method. 

Matrix-element methods are normally based upon two 
stages. First, the detector level kinematics of the events 
of interest are translated into parton-level kinematics via 
transfer functions which describe the effects of the frag¬ 
mentation, hadronization and detector reconstruction. Sec¬ 
ond, the probabilities of predicting a specific parton-level 
event configuration when assuming one or more theoret¬ 
ical hypotheses are calculated on the basis of the associ¬ 
ated matrix elements [32H11I • The simplified method used 
in this work relies on the fact that the event kinematics 
are entirely defined by the configuration of the four top 
quarks, so that one could compute the probabilities using 
this information only. The decays of the top quarks and 
subsequent fragmentation, hadronization and reconstruc¬ 
tion of the final state objects are identical for the sgluon- 
induced and the SM contributions, so that these effects are 
not simulated in detail. Instead, the four-momenta of the 


four top quarks are smeared by sampling a Gaussian distri¬ 
bution with a width corresponding to 20% of the original 
values. As the events are uniquely defined by single points 
in phase space, the CPU-intensive phase-space integration 
normally required by matrix-element methods is avoided. 
Moreover, our method is further simplihed by evaluating 
the probabilities in the SM hypothesis only. 

In order to probe how SM-like are titi events induced 
by the decay of a sgluon pair, we make use of the Mad- 
GRAPH5_aMC@NLO program to generate one sample of 
SM titi events and a collection of samples including sgluon 
contributions for different mg values in the [350, 900] GeV 
range. We then estimate the probabilities in the SM hy¬ 
pothesis of the events in each sample, and compute a quan¬ 
tity known as the discrimination signihcance {D) [45] that 
is used to assess the similarity of the probability distribu¬ 
tions for SM and sgluon events. It is defined as 

jj _ I-PSM ~ -PSgluonl 

where Psm and Psgiuon are the means of the probability 
distributions for the SM and sgluon cases, respectively, 
and where and Pj^^on root-mean-squares of 

these probability distributions. The numerator consists of 
the difference between the means of the two distributions 
in question, while the denominator is the effective resolu¬ 
tion a in measuring this difference. A D value of unity 
would hence correspond to means differing by Icr. 

We have found that in the case of the generated sgluon 
samples, the discrimination significance has a minimal value 
of P « 0.1 that is reached for ms = 400 GeV. The D quan¬ 
tity then rises approximately linearly to a value of P « 0.5 
for ms = 900 GeV. Furthermore, the dependence of P on 
at is found to be negligible. Consequently, the titi kine¬ 
matics for the SM and in the presence of top-philic sgluons 
lighter than about 1 TeV only mildly differ, so that the pa¬ 
rameter space regions of our model in which the titi cross 
section is predicted larger than 32 fb are excluded. 

3.3. Interpretation 

We present in Figure the regions of the top-philic 
sgluon parameter space that are excluded by the dilepton 
and single lepton analyses of Section |3.1| and Section |3.2[ 
The results (solid lines) are represented in the {ms, at) 
plane for a fixed value of Og/A = 1.5 x 10“® GeV“^ (see 
Eq. 0). We then vary Og by 10% up and down, and show 
the induced variations on the constraints by dashed lines. 
For a sgluon mass of about 400-550 GeV, values of at down 
to about 0.6 X 10“^ are excluded. This is related to the 
signal cross section a{pp —> PP —> titi) which is maximal 
in this parameter space region. For smaller and larger 
values of ms, the sensitivity worsens due to the decreasing 
sgluon branching ratio to a top-antitop system and the 
decreasing inclusive sgluon pair cross section, respectively. 

The dilepton analysis turns out to be more constraining 
than the single lepton one, with constraints on ms ranging 
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Figure 2: The excluded regions in (ms,at) space derived from the 
single lepton analysis (red solid line) and dilepton analysis (blue solid 
line). In both cases, the dashed lines correspond to the exclusion 
regions obtained when ag is varied by ±10%. 


ag/A = 1.5 X 10“® GeV“^, a value again typical of su¬ 
persymmetric setups with superpartners around the TeV 
scale. 

Acknowledgments 

This work has been supported in part by the Belgian 
Federal Science Policy Office through the Interuniversity 
Attraction Pole P7/37, by the Strategic Research Program 
“High Energy Physics” and the Research Council of the 
Vrije Universiteit Brussel. Further support has been pro¬ 
vided by the Fonds Wetenschappelijk Onderzoek through 
the “Odysseus” and “Pegasus Marie Curie Fellowship” 
programmes, by the European Commission funding chan¬ 
nel “Marie Curie Intra-European Fellowship for Career 
Development”, by the Science and Technology Facilities 
Council in the United Kingdom and by the Theorie-LHC 
France initiative of the CNRS/IN2P3. 


up to about 750 CeV for a* « 2x 10“^. This analysis has in 
addition the advantage of being not sensitive to the signal 
kinematics, in contrast to the single lepton one. We how¬ 
ever recall that for the parameter space regions that are 
reachable with LHC collisions at y/s = 8 TeV, the kine¬ 
matics differences between the SM and the sgluon cases 
cannot be pinned down. Nevertheless, the situation could 
be different for the LHC Run H expected to probe much 
larger sgluon masses, as the discrimination significance D 
rises with ms (see Section 3.2). 


4. Conclusion 

Many new physics theories predict the existence of new 
colored scalar fields, named sgluons, lying in the octet rep¬ 
resentation of QCD. These fields are expected, at least 
in top-philic sgluon scenarios, to either decay to a pair 
of gluons or to a top-antitop system. In this work, we 
have adopted a simplified model approach describing the 
dynamics of such top-philic sgluons and probed the asso¬ 
ciated parameter space (spanned by one mass parameter 
ms and two coupling strengths a* and ag/A) by extracting 
constraints from two recent analyses of tttt events by the 
CMS Collaboration in the entire LHC collision dataset at a 
center-of-mass energy of 8 TeV. We have found that sgluon 
masses ranging up to 750 CeV are excluded for sgluon cou¬ 
plings to the top quark of a* = 1.5 x 10“^, a value typical 
of non-minimal supersymmetric scenarios with superpart¬ 
ner masses of about 1 TeV or 2 TeV. In the case of smaller 
top-sgluon couplings of at = 0.75 x 10“^ (the smallest 
value which the LHC Run I is sensitive to), the sgluon 
mass is imposed to satisfy ms ^ [400, 550] CeV. We have 
finally also observed that the picture does not drastically 
change with respect to 10% variations of the sgluon-gluon 
coupling strength whose reference value has been fixed to 
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